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ABSTRACT 


A guided missile is one which is usually fired in a direction approximately to the target 
and subsequently receives steering command from the guidance system to improve it’s 
accuracy .The missile being designed has both unguided as well as guided phase. During 
unguided phase, the missile is designed to achieve level flight at a height of 500 m up to a 
range of 3000 m. The guidance system becomes operational to engage the target at 
around 5000 m .In this work, tail deflection required to satisfy the trajectory during 
unguided phase has been postulated using aerodynamic coefficient generation through 
theoretical and wind tunnel methods. The guidance law based on proportional navigation 
has been used in this study. Mathematical model using non-linear aerodynamics has been 
proposed to calculate angle of attack and acceleration response. In order to reduce 
airframe response time, suitable scheme based on Stability Augmentation System (SAS) 
has been proposed. 
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Chapter 1 


INTRODUCTION 

A missile in a simple definition may be called as a self-propelled, unmanned 
space or air vehicle canying an explosive warhead. A guided missile’s path can be 
adjusted during flight, either by automatic self-contained controls or remote human 
control. Guided missiles are of various types and ranges. Missiles may be aerodynamic, 
i.e. controlled by aerodynamic surfaces and following a straight-line trajectory to the 
target, or ballistic, i.e. powered during flight and following a parabolic trajectory. 

Guided missiles are of various types and ranges; long-range missiles generally 
have nuclear warheads, while short-range missiles usually have high-explosive warheads. 
Aerodynamic missiles are of four types. Surface-to-Air and Air-to-Air missiles 
supplement anti-aircraft guns and are often guided by self-contained controls that detect 
and engage the target toward heat or electronic sources. Air-to-Surface missiles, launched 
by aircraft against ground positions, are often radio-controlled. Surface-to-Surface 
missiles (including ship and submarine launched versions) include many different types. 
All long-range missiles are ballistic; the intermediate-range ballistic missile (IRBM) can 
reach targets up to 1,500 nautical miles away, while the intercontinental ballistic missile 
(ICBM) has a range of many thousands of miles. 

Guided missiles apparently had their origin in Germany. For example the Hs.298 
was one of the series of German Air to Air guided missiles developed by Henschel 
Company during world war II l . The Hs.298, which was radio controlled from parent 
aircraft, was to be released either slightly above or below the target. Apparently the 
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height differential made it easier to aim and guide the missile. On December 22, 1944 
three missiles were test flown from JU 88 G aircraft. All three tests were a failure. The 
Rheintocter (R-l) was a Surface to Air missile, developed in Germany during World War 
II. This unusual looking two stage radio controlled missile weighed nearly 4000 lb and 
had three sets of plywood fins: one for the booster and two for the sustainer. The missile 
was ineffective because allied bombers, which were the R-l’s intended target, flew above 
the range (20,000 ft) of this surface to air missile. 

Although the Germans apparently knew the proportional navigation during the 
World War II, no applications on the Hs .298 or R-l missiles using proportional 
navigation were reported 2 .The Lark missile, which had its first successful flight test in 
December 1950, was the first missile to use proportional navigation. Since that time 
proportional navigation guidance has been used in virtually all the world’s tactical radar, 
Infra Red (IR) and television (TV) guided missiles 3 . The popularity of this interceptor 
guidance law is based on its simplicity, effectiveness and ease of implementation. 
Apparently C.Yuan and others first studied proportional navigation at RCA laboratories 
during World War II under the auspices of U.S Navy 4 . 

Guidance is the process of guiding the path of an object towards a given point, 
which may generally be moving. If the given point, which we call as target, is fixed, e.g., 
a seaport then the process is usually called navigation. If the target moves in a way that is 
not predictable -for example, an aircraft evading ground to air missiles-then the process 
is Guidance in it’s narrower sense. 

The guidance law was conceived from physical reasoning and equipment 
available at that time. Proportional navigation was extensively studied at Huges Aircraft 
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Company 5 and implemented in a tactical missile using pulse radar system. Finally a 
proportional navigation was fully developed at Raytheon and implemented in a tactical 
continuous wave radar homing missile 6 . 

After World War II, the U.S work on proportional navigation was declassified and 
first appeared in the journal of applied physics 7 . Mathematical derivation of the 
‘optimality’ of the proportional navigation came more than twenty years later . 

Theoretically the proportional navigation guidance law issues acceleration 
commands, perpendicular to instantaneous missile-target line of sight, which are 
proportional to the line of sight rate and closing velocity. In tactical radar homing 
missiles using proportional navigation guidance, the seeker provides the line of sight rate 
information. In tactical missile proportional navigation guidance are usually implemented 
by moving the fins or other control surfaces to obtain the required acceleration. 

Conventionally, the artillery tank gun uses unguided projectile to neutralize 
enemy tank. The operational range for such tank is well with in 3 km. The projectile or 
shell fired from the tank gun is subjected to higher level of acceleration, non-standard 
atmospheric conditions and thus many times the accuracy in terms of ‘first shot kill’ 
deteriorates. 

From last ten years, considerable efforts are being directed to develop guided 
anti-tank ammunition. Due to strategic reasons no open literature is made available for 
ready reference. 

Armament Research & Development Establishment (ARDE), Pune has 
undertaken a task of developing anti-tank guided missile. The subject missile is being 
developed to meet the following requirements: 
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1. Given maximum elevation of around 10 deg, the missile should be able to reach a 
height of around 500 m within a minimum possible time. 

2. At around 500 m, the missile should maintain almost level flight up to 3 km of range. 

3. At around 3 km range, once the guidance system is made active, the missile should be 
capable of engaging a target at a maximum distance of 2 km using appropriate guidance 
control law technique. 

The first requirement demands large turning rate of the missile. Since during this 
phase booster will be on, part of this turning rate will be obtained from thrust and part 
from lift force generated by the missile at a particular angle of attack. To utilize part of 
the thrust and lift force, it is necessary to introduce angle of attack to the missile. The 
control fins at the rear end of the missile can be appropriately deflected to generate the 
required angle of attack. Since during this phase no guidance and control are operative, 
tail deflection has to be pre-programmed as a function of time to generate required angle 
of attack. 

The second requirement to maintain level flight at around 500 m demands that the 
missile must maintain a predefined angle of attack to balance the weight up to a range of 
3 km. Here again the tail control deflection has to be pre-programmed to generate 
adequate angle of attack to produce desired lift. 

To meet the third requirement, the missiles are configured in such a manner, that 
it generates sufficient acceleration per unit tail deflection to steer the missile towards the 
target. 

The strategic missile with terminal guidance needs to have a high level of 
maneuverability capability at terminal end. In aerodynamic sense it means that the 
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missile should have marginal static stability at the operational angle of attack. The missile 
in question does not have an attached wing as lifting surface and during terminal phase 
thrust is not available, thus most of the lift is to be generated by the blunt shaped 
cylindrical body. The cylindrical body at high angle of attack, generate complex vortex 
pattern to add non-linearity to the flow 9 . The fin attached to the rear end of the missile 
frequently falls in the ‘vortex sheet’ created by the body. Interaction of these vortices 
with the flow near the fin alters the lifting characteristics significantly. More importantly 
physical phenomenon governing this interaction is not well understood or modeled 9 . 
Further at high angle of attack, the body is expected to shed asymmetric vortices, which 
drastically alters the flow field around the fin. Available theoretical methods find it 
difficult to predict this behavior and thus wind tunnel testing remains the best source to 
capture the flow non- linearities and its effect in forces and moments experienced by the 
missile 10 

Accurate values of the forces and moments are a paramount importance to pre- 
program the tail deflection to achieve the desired trajectory and evolve efficient control 
law to implement guidance command. Of course need for an approximate theoretical 
estimation of these forces and moments is absolutely necessary to freeze initial 
parameters of the missile. To refine the design parameters, wind tunnel testing is 
routinely done. Thus as a part of this thesis, aerodynamic parameters were estimated 
using the both wind tunnel testing and theoretical methods. 

The full scale model fabricated by ARDE, Pune was tested in National Wind 
Tunnel Facility at IIT Kanpur.To ensure, sufficient rigidity to the sting, supports mounts 
were fabricated as a part of this exercise. The wind tunnel test data was generated for 
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different combinations of angle of attack and tail deflection. Finally comparisons were 
made between theoretical and wind tunnel estimates. It was generally observed that 
theoretical and wind tunnel estimates compares fairly well as long as angle of attack was 
kept between -15 deg to +15 deg. The wind tunnel data clearly brought out the effect of 
shedding vortices at high angle of attack, which as expected theoretical methods failed to 
capture. 

Once the aerodynamic parameters were frozen based on this study, next step was 
to formulate the trajectory model. Trajectory model was formulated using six degree of 
freedom model to predict range and height of the missile for a given elevator deflection 
and initial conditions. Different tail (elevator) deflections were given as input to the 
trajectory model to obtain the desired trajectory. 

As discussed earlier, to implement proportional navigation, kinematics of the 
problem demands specified amount of accelerations to be made available. Thus, at the 
design stage, the normal acceleration produced per unit tail deflection is obtained using 
short period approximation 11 . Non-linear model 10 was assumed to calculate steady state 
angle of attack and normal acceleration. Next, the normal acceleration obtained was 
compared with acceleration required as per proportional navigation kinematics to freeze 
the amount of tail deflection required for terminal guidance. The normal acceleration is 
computed assuming the steady state value of a i.e. (a mm ) for particular tail deflection 5 . 

However since these missiles are designed to have low static stability (for high 
maneuverability), the pitch damping is very low and thus takes longer time to reach 
steady state value of a tnm to generate required normal acceleration. This is highly 
undesirable from implementation point of view. Thus a need is arised to artificially damp 
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the missile dynamics without altering the steady state trim value. Hence, as the final work 
of this progress, preliminary stability augmentation scheme has been proposed for 
implementation. 

Chapter 2 discusses in detail the theoretical methods employed and the conduct of 
wind tunnel testing. Chapter 3 discusses strategy development to pre-program tail 
deflection to satisfy the first missile requirement. Chapter 4 presents scheme to engage 
target using proportional navigation. Further, theory employed to generate accelerations 
through tail deflections has been discussed in detail. It also discusses the method to 
design artificial stability augmentation system (SAS) to artificially damp the missile 
dynamics. Finally, in chapter 5 the conclusion and scope for fixture work have been 
discussed. 
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Chapter 2 


ESTIMATION OF AERODYNAMIC PARAMETERS 
2.1 GENERAL 

The conventional approach hitherto for understanding the flight behavior of 
missiles is to develop mathematical models that could predict all elements of the 
trajectory from launch to target. To this purpose it becomes essential that all the forces, 
moments affecting the flight of the missile be accounted in a well-defined mathematical 
form 11 . Beginning with the most simple but relatively inaccurate model, the in- vacuo 
trajectory model, more and more sophisticated models of increasing accuracy such as 
point mass model (PMM), the modified point mass model (MPMM) and six degrees of 
freedom model (SDFM) have been developed. These models require aerodynamic 
coefficients as input. 

There are three distinct approaches for estimating aerodynamic parameters: 

1 . Theoretical methods 

2. Wind Tunnel Testing 

3. Flight Testing 

At the primary design stage of any system, theoretical methods 10, 12,13 are useful in 

spite of their limited accuracy. The wind tunnel testing improves the accuracy of 

estimation but it is time consuming and expensive way of estimating the aerodynamic 

parameters. Precise simulation of control surfaces, power effects, and flight condition is 

difficult. The model tested in the wind tunnel is generally slightly different from actual 
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flight due to last minute configuration changes. Other reasons of discrepancies between 
flight and wind tunnel results are Reynolds number discrepancies and interference due to 
support system. It is therefore, desirable that the wind tunnel estimates be corroborated 
with the estimates from actual flight test data. 

In the present work, the aerodynamic coefficients were obtained using theoretical 
methods and wind tunnel testing. The model being tested is having a blunt nose with 
cylindrical after body. The fins attached to the rear part of the body is having tapered 
configuration. There are also noticeable gaps between the fin base and cylindrical 
surface. The missile being a tail controlled (no wing) needs large angle of attack to 
generate sufficient normal lift for maneuver. Hence, in general the forces and moment 
generated on the missile will be highly non-linear in nature. Further, at high angle-of- 
attack, the body will shed vortices and subsequently these vortices will interact non- 
linearly with flow field around the tail 9 . The theoretical method needs to be powerful 
enough to capture both linear and nonlinear behavior. Any theoretical methods will find 
it difficult to estimate aerodynamic parameters accurately for such model configuration. 
Thus, extensive wind tunnel testing is necessary to arrive at the final configuration. 
Details of the theoretical methods and wind tunnel testing that were used to estimate 
force coefficient are given in the subsequent paragraph. 
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2.2 DESCRIPTION OF THEORETICAL METHODS TO ESTIMATE FORCE 
AND MOMENT COEFFICIENT 
2.2.1 BODY ALONE 

There are many non-linearities that occur in the weapon aerodynamics. The ones that 
have most influence on the body alone are angle of attack, mach number, crossflow 

reynolds number and asymmetric vortices. This entire phenomenon can be properly 

/ 

modeled in an approximate sense except asymmetric shedding vortices. 

Normal force coefficient can be written as 


(C„) =(C N ) +(C N ) (2.1) 

V N J Body v N ' Body, lin ' N J Body, nonlm V ' 

(C N ) B m ca n be approximated using munk factor 12 (k 2 -k^). However this approach 

is strictly valid for slender body with lifting nose. 

The missile under investigation has slenderness ratio of around 8.33 with blunt 
nose. The theoretical estimates are going to be inferior for such a “blunt nose slender 
body”. Due to non-availability of any other methods, (C N ) body lm has been computed 
using munk factor obtained through Fig. 2.2. 1 below. 




M<1.0 



FINE NESS RATIO I/D 


Fig. 2.2.1 Variation of Munk Factor with fineness ratio 

10 




(C N ) is a cross flow term based on drag force experienced by an 

Body,nonlm 

element of circular cylinder of same diameter in a stream moving at the cross 
component of the stream velocity V m sin a .The cross flow term is primarily created by 
the viscous effects of the fluid as it flows around the body, often separating and creating a 
nonlinear force coefficient. The non-linear coefficient ( C M ) can be approximated by 


where 


(C„) =rjCd c .-^-. sin 2 a 

N v Body, nonlm ° A j- 




(2.3) 

(2.4) 


A p =Plan form area; A ref =Reference area 


M n = M sin a (2.5) 

t] used in the above expression is the drag proportionality factor or the ratio of crossflow 
drag of a cylinder of finite length to one of infinite length. This variation 10 is graphically 
presented in Fig. 2.2.2. 



4 8 12 16 20 24 28 32 36 40 

Circular Cylinder Lengtft-to-Diameter Ratio 
Flat Plate Lenath-to-Width Ratio 


Fig. 2.2.2 Ratio of cross flow drag coefficient for a finite length cylinder (or flat 
plate) to that for infinite length cylinder (or flat plate) at Mn = 0 
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Procedure for Cd c calculation: 

The cross flow drag coefficient Cd c cm be computed as a function of M N using 
Fig 2.2.3. 



Fig. 2.2.3 Variation of cross flow drag coefficient with cross flow Mach number 


It is important to identify the nature of flow i.e. sub critical or super critical for 
estimation of cross flow drag coefficient Cd c . Using the following conditions could be 

used to identify the nature of the flow. 

Super critical: 


If R„ > R 


'AW NC 




Sub critical: 


If < R nc 


M n <M nc 


where Rnc = 180,000 and Mnc = 0.1 have been assumed 10 . 
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Thus knowing the values of 77 and Cd c for a given mach number and angle of attack, C N 
could easily be calculated using eq. (2.1). 

To locate the point of application of this force, the center of pressure at subsonic 
speed has been assumed to be at the centroid of the planform area 10 ' The center of 
pressure moves aft /forward with angle of attack/mach number. The variation of center 
of pressure for body alone can be computed using the data base given in Ref.( 10 ). The 
database used for this case is presented in tabular form in table 2.2.1 for ready reference. 


Table 2.2.1 Shift in body alone center of pressure as a function of Mach number 
and angle-of-attack (as a fraction of body length) 


Shift in body-alone center of pressure as a function of Mach number and AOA (as a fraction of body length) 
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2.2.2 FIN ALONE NORMAL FORCE COEFFICIENT 


The fourth order equation for the fin alone normal force was found to be most 
accurate for all angle-of-attack 10 .The fin alone normal force coefficient is thus defined as 
0 C N )fm = a o +a x a f +a 2 a 2 f + a 3 a 3 f + a 4 a* f (2.6) 

where 

a f =\a + 6\ and S is the fin setting angle. 

To predict the fin alone normal force coefficient, the equation requires these five 
constants to be evaluated. Thus five independent equations or conditions are needed. The 
first condition to evaluate a 0 is that most weapon lifting surfaces are symmetric and have 

zero camber. As a result normal force at zero angle of attack is zero. 

Thus 


a 0 =0 

a \ ~ (^Na Lo 

which is evaluated using 10 


(2.7) 

(2.8) 


«,=(C„„U= — r (2.9) 

2+[aR 7 (P 1 -Han 1 A,„)+4l 2 

where 

/? = -\/ 1 - M 2 

The other coefficients a 2 , a 3 and a 4 can be evaluated using the following expressions 10 

a 2 = 34.044(C„) ii , 1J , -4.824 (C„)„ JS . + 0.426(0,,)^ -6.412a, (2.10) 

a, =88.240(C M ) a . rf +23.032 (C„)„ JJt -2.322(C„)^ tf +11.4640, (2.11) 

a, = 53.219(C w ) a _ ]j0 -17.595(C„ +2.661(C„) ot<0 . -5.971a, (2.12) 
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The terms (^v) a=15 o >(C^) a=35 „ ,(C w ) a=6o0 are the normal force coefficients at 

cc f = 15°, 35°, 60° respectively as a function of Mach number, taper ratio and aspect ratio. 

The numerical value of these coefficients as function of Mach number, taper ratio and 
aspect ratio are tabulated in Appendix A, Table 2(A). 


2.2.3 FIN BODY AND BODY FIN INTERFERENCE 


There are many types of interference effects that can occur in aerodynamics. First 
of all, when a fin body configuration is present, there is an additional pressure load 
induced on the body as a result of fin being present and likewise, there is an additional 
pressure load induced on the fins as a result of being in the presence of body. Another 
type of interference is that due to the body shed vortices as they affect the fin surface. 

There are two primary types of interference K f(B) and K B(f) . These are 

interference factors associated with normal force of the fin in the presence of the body 
and additional normal force on the body as result of fin being present due to angle of 
attack. 

Mathematically 


K 


'N 


f(B) 


f(B) 


(2.13) 


(^N ) fin 

and defined as Ratio of normal force coefficient of fin in presence of body to that of fin 
alone at 6 = 0 deg. 


AC 


K 


N 


B(f) 


B(f ) 


(Cv) 


(2.14) 


NJfm 
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and is defined as ratio of additional body normal force coefficient in the presence of fin 
to fin alone normal force coefficient at S = 0 deg.The mathematical model to estimate 
K f{B) and K S(f) are presented next for ready reference 10 . 


MATHEMATICAL MODEL FOR K f(B) 


K f(B) “[^/(B)] ssr + 


L=o ; 


*/<»=[* 


[a£/ (B) _ 

ft 

II 

o 

+ 

III 

-a c ). 

Kf(B) = 1 “ 

K a M ~a D K 

fw\=a Dt 

h 

*,<«=[* 

K f<n\ D ~ 


max 


da 

a D 

~ a c 




dK 




da 


a<a c 
a c <a<a D 

a D <a<a M 

a>a M 


[ x /<«L=[ r /<»>L + [ AK /«>L 


MATHEMATICAL MODEL FOR K, 


B(f ) 


K *(f) ~ [^C/)]ssr + [^(/>L + M* 


dK 


f(B) 


da 


a<a. 


K «</) =K/>L, ; “■ <a -■ ** 


k b(/) ~ \_ k bu) ] min ’ a ~ a i 


(2.15) 

(2.16) 

(2.17) 

(2.18) 

(2.19) 

( 2 . 20 ) 

( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 


16 



Jw 


< r" 
1 + - 

V S J 


/ \2 
' r 


— tan -1 — ^ 
2 2 





r l 

f — - — l + 2tan _1 — 

) 4 I 

5 2 

\r sj s_ 


V sj 


A 2 


1 -- 




^SBT ~ + s J \. K m ] 


SBT 


(2.24) 


(2.25) 


The data base to estimate [a K f(B) , a c , [K f(B) ,a D ,a M , [a K B{f) , 


dK 




da 


,a,. 


a 2 for different combination of is a function of aspect ratio, taper ratio and Mach number 
are enclosed in Appendix ‘A’, Table 2(B). 

Thus the normal force coefficient for the complete body and fin incorporating 
interference factors can be written as: 


c N HC N ) Body HC N ) f AK fW+ K B(/) ). 


* fin 

u 


(2.26) 


CENTER OF PRESSURE OF MISSILE 

Total center of pressure of the missile is estimated using the following relation 

dC t + C Nnorilm .X 2 + C N j m X 3 


(Xcp) 


missile 


c +r + c 

^ NJrn ^ ^ N,nonlin ^ ^ N Jin 


(2.27) 


where 


X l ^Center of pressure for nose linear loads «0.6 nose lengths. 

X 2 =Center of pressure of non-linear body load .The non-linear center of pressure shifts 
with angle of attack. 
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X z =Center of pressure for fin, is assumed to be at the quarter chord point(c/4). Its 
variation with angle of attack has been neglected. 

Total moment coefficient has been obtained by taking moment about the missile center of 
gravity. 



[C„ ) + C 


"AT 


.(X 2 -X eg ) + C^.(JT 3 -X^)] 

d 


(2.28) 


2.3 DESCRIPTION OF WIND TUNNEL TESTING TO ESTIMATE FORCE AND 
MOMENT COEFFICIENTS 

In any wind tunnel testing, the accuracy of wind tunnel estimates strongly 
depends on the exactness of the model, precise mounting to simulate test conditions and 
of course accurate calibration of data acquisition system. 

Armament Research & Development Establishment (ARDE), Pune , supplied the 
full-scale model for design evaluation, analysis and exhaustive wind tunnel testing. The 
schematic of the model used for testing is presented in Fig. 2.3.1. The major dimensions 
of the model are not mentioned for the sake of secrecy. 

1 

d 

I/d=8 33 V_7 

r/d=0 5 

Fig 2.3.1 Schematic of the model used for design evaluation and testing. 
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The objective behind the wind tunnel testing of this model was to estimate force 
and moment coefficient at speed of around 60m/s for various combination of angle-of- 
attack (a) and tail setting angle (8). For a fixed tail settings, the angle-of-attack was 
varied from -15° to 45°. Tests were repeated for various angles of tail setting angles 
ranging from +10° to -40°. ‘+’ is the tail down and ‘ - ’ is the tail up movement. 

The full-scale model was tested at National Wind Tunnel Facility (NWTF), at IIT 
Kanpur. The sting available in this tunnel was not sufficient to give required large angle- 
of-attack by rotating the model in vertical plane. Hence, a new methodology was adopted 
for testing the model. Angle-of-attack was simulated using two vertically located 
turntables in the test section. The model was placed between these turn tables. A fitment 
as described in Fig. 2.3.2 was fabricated to hold the model rigidly between these turn 
tables. To avoid interference due to wake formation by the fitment, the model was 
positioned at a sufficiently large distance from the vertical element of the fitment. 

A six-component balance was installed inside the model to measure forces and 
moments. The tunnel was stabilized at a wind speed of 60 m/sec and data acquisition 
system was switched on to acquire data for three missile configuration namely, body- 
alone, body with fins at zero settings in *+’ configuration, and body with fin at different 
positive (down) and up (negative) fin setting angles. Before every run, a dry run (no 
wind) was conducted to estimate the bias error, if any. The raw data acquired were then 
processed to obtain desired force and moment coefficients namely C x ,C Y ,C N ,C m ,C n ,C,. 
The detail procedure followed for processing wind tunnel raw data is given in 
Appendix“B”. 
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Fig 2.3.2 Schematic of fitment fabricated to hold the model. 
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2.4 RESULTS AND DISCUSSION ON ESTIMATION OF AERODYNAMIC 
PARAMETERS USING THEORETICAL AND WIND TUNNEL METHODS 
2.4.1 WIND TUNNEL TESTING 

The force and moments coefficients were obtained using wind tunnel data. The 
wind tunnel raw data was processed, corrected for bias and then corrected to non- 
dimensional form to get aerodynamic coefficients C x ,C Y ,C N ,C m ,C n ,C,. 

Wind tunnel data for three-missile configurations-body alone, body with fins at zero 
settings in *+’ configuration, and body with fin at different positive (down) and negative 
(up) fin settings are analyzed in this section. All these configurations were used to record 
data with a varying from -1 5° to 45°. 

The Normal force ( C v ) and Pitching moment ( C m ) variation with angle-of-attack 
for varying fin angle settings ( S ). 

The Normal Force: For body alone. Fig. 2.4.1 shows that Cm Vs a variation is nearly 
linear in the range of a from -15° to 45°. A positive value of C N is observed at a = 0 
while the corresponding C m is negative (-0.01 15). Because the center of pressure of body 
is ahead of c.g, a positive value of Cm and negative C m are inherently inconsistent. To 
check our conjecture that this inconsistency is due to error in zero setting of a, a plot of 
Cm for a -5° to 5° was plotted and joined by a straight line. The intersection of this 
straight line with Cm axis at a = 0 showed Cn = -0.018. Thus, the Cm and C m values 
are consistent in sign: A negative Cm ahead of c.g producing a negative C m . 
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Fig. 2.4.1 Normal force coefficient for body alone configuration 

It may be mentioned that the a = 0 setting for body plus fin combination was 
redone afresh, it seems that the new setting of a yielded Cm (-0.098) and C m (0.0147) 
values, which are consistent This further suggests that the inconsistency of Cm and C m 
for body alone was probably due to a = 0 setting only. 

Civ Vs a at <?= 0: Figure 2.4.2 shows that the C N Vs a is almost linear till o=10°, and 
then non-linear contribution is almost as much as the linear contribution till «=20°. At 
o=20°, there seems to be break in the trend and a new non-linear behavior is observed. 
At a=10°, 5=0 Cm was measured to be 0.829. Also, the measured value of Cm for body 
alone at «=10° was 0.555. These numerical values will be compared next with the 
measured values at different a and 5 combinations. 
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Fig. 2.4.2 Normal force coefficient for body and fin at cant = 0 deg. 

Ov Vs a at 8 ± 0: Figure 2.4.3 shows the trend for 8 = -10° . It has a trend of CV 
variation with a similar to the one observed for 8 = 0. 



Fig. 2.4.3 Normal force coefficient for body and fin at cant = -10 deg. 

The measured Cat at o=10° and 8— -10° was 0.797. One would expect that at this 
combination of a and 8 ,fin contribution to is zero and Cat would have value equal 
to that for body alone at o=10° i.e. Cn =0.555 .However , the measured value is 
actually comparable to 0.829 measured for body plus fin at a= 10° .It is known that for 
long slender bodies, flow separation and shedding of trailing vortices begins from the 
body at some distance from the nose of the body and this vortex flow would interact with 
fins to generate normal force .This we believe explains the observed value of Cat The 
variations of measured C.v with a for 8= -20° ,-25° ,-30° are presented in 
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Fig.2.4.4-2.4.6. Again, as in case of a =10° and 8 - -10, the combination of (a=20° and 
5- -20°), ( a =25° and 8= -25°) and ( a =30° and 8— -30°) show that measured values 
of Cat values greater than that for body alone at the corresponding a. Here again, it seems 
to confirm our earlier conjecture that vortex flow is responsible for additional Cm -In an 
overall way, the variation of Cm with a for all these 8 settings are similar in character as 
discussed for 8 =0 case. 



Fig. 2.4.4 Normal force coefficient for body and fin at cant = -20 deg. 

a, deg 



Fig. 2.4.5 Normal force coefficient for body and fin at cant = -25 deg. 
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Fig. 2.4.6 Normal force coefficient for body and fin at cant = -30 deg. 

The Pitching moment: 

For body alone. Fig 2.4.7 shows that the variation of C m Ys a is increasing almost 
linearly up to a = 10°, and then non-linearly till a- 30°, beyond which the C m decreases. 



Fig. 2.4.7 Pitching Moment coefficient for body alone configuration. 

It was noted earlier that the normal force (body) increases with a even beyond 
a=30°. This would suggest that the decrease in C m is due to a shift of the center of 
pressure closer to C.G. 

Cm Vs a for S=0: Figure 2.4.8 shows that the missile is stable till a=T0° but as a varies 
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between 10° to 30°, the missile attains almost neural stability and then again beyond 
a=30°, it becomes stable (more stable than at o<10°). 



Fig. 2.4.8 Pitching Moment coefficient for body and fin at cant = 0 deg. 

It was pointed out that body alone contribution increases till a =30° and then 
drops off. The reason for achieving neutral stability in the range of 10° < a < 30° seems 
to be a result of well designed tail size and its location. For efficient maneuvering, it is 
desirable that the missile has a low or almost neutral stability in the range of a that is 
likely to be used during maneuvers. The C m contribution from the fin and body are made 
to cancel each other; the C m (negative contribution) due to fin increases as a increases but 
so also the body contribution (positive) increases with a. Now by clever design of fin, it 
is ensured that both contributions cancel each other. This is achieved for values of a up 
to 30°. However, as mentioned earlier for a > 30°, the unstable contribution to C m due to 
body drops off and thereby the net C m again is increasingly negative as a increases, and 
missile becomes stable again. 
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C« Vs a for 5 _ < 0°: . 



Fig. 2.4.9 Pitching Moment coefficient for body and fin at cant = -10 deg. 

Figure 2.4.9 shows C m Vs a for S= -10° .It is noted from the figure that missile is 
neutrally stable in the range of 15° < a < 30°. For a = 10° 8= -10°, one would expect to 
measure C,„ equal to values of body alone, if the fin experienced no normal force. But as 
pointed out earlier, the shed vortex flow from body seems to give rise to normal force on 
fins and this in turn, will create a nose down C m .As a increases ,the vortex flow at the 
corresponding a would pass over the fin increasing distance and less and less normal 
force would be induced. However, increase in a will increase the fin contribution of Cm 
due to increase in angle of attack seen by the fin. These two effects seem to add to 
increase C N with a in such a way that the negative C m produced is more in magnitude as 
compared to the positive C m due to body in the a values varying between 15° to 30°. The 
difference in body and fin contributions remains a constant in the range of 15° < a< 30° 
and thus the observed neutral stability in this range of a. Also it is noted that at o=10°,the 
negative C m due to fin and positive C m due to body exactly cancel each other, providing 
the condition (C m =0) for a = 10° and 5 = -10° .Beyond 30°, again due to 
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drop off in body contribution, the overall C m becomes increasingly negative with a and 
thus makes the missile stable. 



Fig. 2.4.10 Pitching Moment coefficient for body and fin at cant = -20deg. 

For 8 = -20°, -25° and -30°, the neutral stability is again observed but not so 
emphatically. For example, Fig. 2.4.10 for 8= -20° shows that C m in the range of 
20°<a<30° is changing from negative at a = 20° to positive at a - 25° and to negative 
again for oc= 30°. If fluctuations were assigned to experimental errors, and C m value is 
averaged out for 20° < a < 30°, the C m Vs a would again show neutral stability. Beyond, 
o=30°, as explained earlier, net C m starts to build up to more and more negative values 
and missile is stable again. Also, as observed for a =10° ,8= -10° combination ,the 
a= 20° ,8= -20° combination also shows trim condition (C m =0) . 
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Body Axes Moment Coefficients 



Fig. 2.4.11 Pitching Moment coefficient for body and fin at cant = -25deg. 

A less than obvious similar trend can be observed by closely studying Fig. 2.4.11 
that shows C m Vs a for 8= -25°. Here again, for 25° < a < 30°, the C m Vs a curve seems 
to relatively flatten out before again becoming negative for a >30°.Also it is noted that 
for a = 25°, 8 = -25° combination, the trim condition (C m =0) is achieved. 

Finally, Fig. 2.4.12 shows results of C m Vs a for 8= -30°. 



Fig. 2.4.12 Pitching Moment coefficient for body and fin at cant = -30deg. 

Here one would expect no neutral stability as we are already at 8- -30° and for e£>30°, 
the decreasing value of body C m would not get cancelled by fin contribution and thus a 
negative C ma is observed. Also, it is noted that unlike for other combinations of a = x° 
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,8= -x° for a = 30° ,8= -30° no trim (C m = 0) is achieved , and a net small value of C m 
is observed, and trim is observed at slightly high value of a - 32° 

Trim of the Missile: As pointed out, it is worth emphasizing that the missile trim (C OT =0) 
is achieved at values a =-S for all the settings of 8 Fig. 2.4.7- 2.4.12 discussed earlier 
show this for all the 8 settings. One case of positive 8 settings of 10° is given in Fig 
2.4.13 Here again, it can be seem that the missile trims at a = -10°. This is a very 
important aspect of the design. The a trm is achievable by simply deflecting the fins to the 
corresponding value of opposite sign. 

Body Axes Moment Coefficients 



Fig. 2.4.13 Pitching Moment coefficient for body and fin at cant = +10deg 
Conclusion: It appears to be a well-designed missile from operations point of view. It 
has neutral or low stability in the desired range of trim angle-of-attack range. For further 
understanding and appreciation of aerodynamic flow past the missile at varying a and 8, 
flow visualization test are planned in near future. Also to scan the data for varying a, a 
new model is being fabricated to meet stringent accuracy requirements and also enable 
the fin settings to be changed by remote control during the tests. The present model 
seemed to have slight curved fins and slight asymmetry in geometric shape. 
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FORCE COEFFICIENT C y 

For body alone, referring to Fig 2.4.14 it can be seen that C r is close to zero for 
a up to 20°. Since the model is symmetric and is not subjected to any sideslip angle, one 
should expect C Y = 0 for all range of a under normal circumstances. However it can be 

seen from the Fig. 2.4.14, that C Y abruptly becomes positive between a = 20° to a = 40° . 
A positive C Y should result in negative C n . The wind tunnel data indeed show negative 
C n for the range of a = 20° to a = 40° . This is attributed to generation of asymmetric 
vortices at high angle of attack. 



a,d*g 

Fig 2.4.14 C Y ,C n from Wind Tunnel as a function of angle of attack for body alone 
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2.4.2 THEORETICAL METHOD: 


Theoretical methods explained in earlier section was applied to estimate force 
coefficient C N and moment coefficient C m . For body alone and C m were estimated for 

different values of angle-of-attack [-15° to 45°]. The results thus obtained were 
compared. 

2.4.3 COMPARISION OF THEORETICAL ESTIMATES WITH WIND TUNNEL 
DATA 

BODY ALONE 

To start with, normal force coefficient C N for body alone was computed using eq (2.2) 
and eq (2.3). Theoretical estimates of C N along with corresponding wind tunnel 
estimates given in Fig. 2.4.15 and Table 2.4.1. 

Table 2.4.1 Comparison of Normal Force coefficient from wind tunnel and 
theoretical estimates for body alone configuration 


Alpha, Deg, 

Theoretical 

Wind Tunnel 

-15 

-1.0261 

-0.97323 

-10 

-0.56284 

-0.51578 

-5 

-0.21944 

-0.25778 

0 

0 

0.12023 

5 

0.21944 

0.22187 

10 

0.56284 

0.55578 

15 

1.0261 

0.96654 

20 

1 6015 

1.3659 

25 

2.2848 

1.6143 

30 

3 0483 

1.8262 

35 

2.5464 

2.2399 

40 

3.0184 

2.7712 

45 

3.4924 

3 028 
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a, deg 

Fig 2.4.15 Comparison of Normal Force Coefficient from wind tunnel and 
theoretical estimates for body alone configuration. 

Column (2) and Column (3) of table 2.4.1 lists the numerical values of body alone normal 
force coefficient C N obtained through the theoretical and wind tunnel methods. It can be 

easily observed that within the a range of -15 to +15 deg, the estimated values of the 
coefficient compare well with the wind tunnel estimates. Beyond 20 deg., difference 
between the C N estimated by the theoretical and wind tunnel method widens 

significantly, especially at angle of attack of 30 deg. This was expected as the theoretical 
methods (available /used) is not powerful enough to capture the effect due to asymmetric 
vortices at high angle of attack. Despite these limitations, in general theoretical prediction 
for the body alone is fairly good and can be used for initial design and analysis of the 
missile dynamics. 
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BODY AND FIN CONFIGURATION 


The fourth order equation for the fin alone normal force as given in equation (2.6) 
was used to compute C N due to fin alone. In applying the fourth order equation, it was 

observed that addition of third form onwards yielded unrealistic value for C N fm Since, 

angle of attack of interest was primarily up to 20-30 deg, it was decided to estimates 
C N fm using first three terms of eq (2.6). Table (2.4.2) to (2.4.7) present tabular 

comparison between estimated value of C N obtained through wind tunnel testing and 
theoretical method for different combination of angle of attack and tail setting angle. The 
values of C v for different cant angles (-10 to -40) and angle of attack (-15 to 40) are 
compared graphically and presented through Fig (2.4.16) to (2.4.22). 
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Table 2.4.2 Comparison of wind tunnel and theoretical estimates for body 
and fin at cant=0° configuration 


Alpha, Deg, 

Theoretical 

Wind Tunnel 

-15 

-1 5258 

-1.2869 

-10 

-0.94077 

-0.8899 

-5 

-0.4308 

-0.49837 

0 

0 

-0 098257 

5 

0.4308 

0.40483 

10 

0 94077 

0.82968 


1.5181 


2.1564 


25 

2.8571 

2.1246 

30 

3.5927 

2.4409 

35 

3.0105 

3 204 

40 

3.3718 

3.8027 

45 

3.7118 

4.3538 


Theoretical 
Wind Tunnel 


— 

Body 

, 

+ Fin Configuration 

— 

_ * 





Fig 2.4.16 Comparison of wind tunnel and theoretical estimates for body 
and fin at cant=0° configuration. 





























Table 2.43 Comparison of wind tunnel and theoretical estimates for body 
and fin at cant=10 0 configuration 


Alpha, Deg, 

Theoretical 

Wind Tunnel 

-15 

-1.4397 

-1.3803 

-10 

-0.56284 

-0.80002 

-5 

-0 21019 

-0.12293 

0 

0 19722 

0.22698 

5 

0.55982 

0.56841 

10 

1.0016 

1.0792 

15 

1 5128 

1.4198 

20 

2.0875 

1.7351 

25 

2 7267 

2.1161 

30 

3.404 

2.5397 

35 

2 7736 

2.923 

40 

3 0914 

4.0674 

45 

3.3923 

4.13 


















Table 2.4.3 Comparison of wind tunnel and theoretical estimates for body 
and fin at cant=-10° configuration 


Alpha, Deg, 

-15 

-10 

-5 

0 

5 

10 

15 

20 

25 

30 

35 

40 

45 


Theoretical 

- 1.5184 

- 1.0016 

- 0.55982 

- 0.19722 

0.21019 

0.56284 

1.4299 

2.1316 

2.8939 

3.6879 

3.1539 

3.5588 

3.9376 


Wind Tunnel 

• 1.3801 

- 0.97491 

- 0.51546 

- 0.28287 

0.13138 

0.79674 

1.2271 

1,8027 

2.0677 

2.3512 

3.0214 

3.6635 

4.1373 






















Table 2.4.4 Comparison of wind tunnel and theoretical estimates for body 
and fin at cant=-20° configuration 



























Table 2.4.5 Comparison of wind tunnel and theoretical estimates for body 
and fin at cant=-25 0 configuration. 
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and fin at cant=-30 0 configuration. 



Fig 2.4.21 Comparison of wind 
fin at cant=-30 0 configuration. 

















Table 2.4.7 Comparison of wind tunnel and theoretical estimates for body 
and fin at cant=-35 0 configuration. 


Alpha, Deg, 

Theoretical 

Wind Tunnel 

-15 

-1.0906 

-1.4313 

-10 

-0.74433 

-1.0391 

-5 

-0.47308 

-0.8554 

0 

-0.28099 

-0.4415 

5 

-0 044081 

-0 18083 

10 

0.36162 

0.04584 

15 

0.92535 

0.6772 

20 

1 6342 

0.90318 

25 

2.4828 

1.423 

30 

3.4346 

1.9322 

35 

2.5464 

2.3987 

40 

3.6168 

3.1598 

45 

4 0931 

3 5773 



Fig 2.4.22 Comparison of wind tunnel and theoretical estimates for body and 
fin cant=-35 0 configuration. 
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A closer look through these tables reveals that up to 15 deg, consistent matching 
exist with the wind tunnel estimates. In fact for =0 the matching between the both the 
estimates are fairly good; However the differences between theoretical and wind tunnel 
estimates widens up as tail angle is increased. In fact beyond angle of attack of 20 deg, 
where in fact asymmetric vortices have highly non-linear influence on fin performance, 
the prediction is poor. Even for small angle of attack as cant angle is increased the 
theoretical prediction shows large noticeable scatter from reference wind tunnel 
estimates. 

Further it was also seen that for a = -8 , the theoretical predictions fail to capture 
the correct values of C N and then C m also. This is well understandable, as the method 
used here cannot take into account of effect of symmetric shedding vortices generated by 
the body on fin in estimation of normal force. 

PITCHING MOMENT COMPARISION 

Figure 2.4.23, Fig 2.4.24 graphically compares the values of pitching moment 
coefficients as obtained from theoretical methods and wind tunnel testing for a body 
alone configuration an d for body and fin at a cant angle setting of zero deg. The 
comparison of C m for body alone is presented in Fig 2.4.23 .It can be observed that 

theoretically estimated C m does not compare well with the measured C m . This is 
attributed to incorrect estimation of center of pressure by theoretical methods. Similar 
trend was observed while comparing estimated and measured C m for body and fin at zero 
setting angle. This is pictorially presented in Fig. 2.4.24.Therefore, the center of pressure 
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V 


as obtained from the wind tunnel estimates was used for calculating the theoretical C m 
and it showed a close proximity with the measured C m as shown in the figure; designated 
with * mark. 



Fig 2.4.23 Comparison of Theoretical Estimates with Wind Tunnel for body 
alone. 



Fig 2.4.24 Comparison of Theoretical Estimates with Wind Tunnel for body 


and fin at cant=0 configuration. 

* - Center of pressure obtained from wind tunnel used. 
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Chapter 3 

TRAJECTORY MODELLING AND DEVELOPMENT OF CONTROL INPUT 
STRATEGY 

3.1 GENERAL 

The six degrees of freedom model is the ultimate in trajectory models allowing 
missile to have yaw in two planes as well as three spatial degrees of freedom and spin. It 
is not necessary to make any assumption regarding the linearised aerodynamics and the 
results obtained from the model would be exact. However the in determinability of many 
of the initial conditions and aerodynamic coefficients which are required as input 
frequently results in model not giving significantly good results. Nevertheless it is 
powerful tool for an ammunition designer. 

3.2 DEVELOPMENT OF SEX DEGREES OF FREEDOM TRAJECTORY 
MODELLING 

The six degrees of freedom model is based on the basic Newton laws which state that 
summation of all external forces in any direction must be equal to the time rate of change 
of momentum and the summation of all the moments of external forces must be equal to 
the time rate of change of angular momentum, all measured with respect to axis fixed in 
space. This can be stated mathematically as 


Yu F x= m - a x 

( 3 - 1 ) 

Yj F y = m.a y 

( 3 . 2 ) 

Yu F ^mxi z 

( 3 . 3 ) 
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(3.4) 


y i _ dH_x 

^ dt 




Z N - 


dt 

dH z 

dt 


(3-5) 

(3.6) 


The angular momentum vector can be represented as 


H = 


"4 I xy a 

~ 4y lyy — 4 

\~^xz ~^yz I a J 


XO 


(3.7) 


The components of co are p,q,r in the body axes system by definition and thus can be 
written as 


co = pi+ qj + rk 

With the above substitutions the forces and moment equations can be written as 
F x =m(u + qw-rv) 


(3.8) 


(3.9) 


F = m(v + ru- pw) 


F z = m(w+ pv-qu ) 


M x = pi XX - ql xy - rl„ + qr(I B - /^ ) + (r 2 - tf 2 )/* “ pql„ + rpl xy 

M y = -pixy + qlyy ~ *1 „ + /J*/* - IJ + ~ T* )I„ ~ qtl xy + P ql ^ 


(3.10) 

(3.11) 

(3.12) 

(3.13) 


K = -pi XX + qlyx + fl a + pq(I„ - /«) + (q 2 - P 2 )lxy ~ r plyx + „ 


(3.14) 


3.2.1 EXTERNAL FORCES 


The external forces F x ,F y ,F z can be written as 


P = p , p . p 

* x ”x,aero ' * x,grav ' " x f Thntsi 


(3.15) 
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F =F +F 

A y y,aero y,grav 

(3.16) 

F : = F :„ + F : x „, 

(3.17) 


(3.18) 

II 

(3.19) 

Kacro^SC z 

(3.20) 

F XRmv = ~mg s'm 0 

(3.21) 

= mg cos <9 sin ^ 

(3.22) 

F :Rrav = mg cosdcos</> 

(3.23) 

The external force equation can now be collectively written as 

F x = qSC\ - mg s\n0 + T 

(3.24) 

F y - qSC v + mg cos 0 sin (j) 

(3.25) 

F, = qSC, + mg cos 6 cos (j) 

(3.26) 

3.2.2 EXTERNAL MOMENTS 

The external moments M r ,M v ,M r can be written as 

M x =qSbC, 

(3.27) 

M y = qSbC m 

(3.28) 

M, = qSbC n 

(3.29) 


It should be noted that the external forces and moments are w.r.t the inertial frame 

of reference i.e. Earth fix axes while the momentum and angular momentum terms are 

w.r.t body axes system. Thus suitable transformations of body fixed axes system to earth 
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fix axes system is necessary as the Newton mechanics is valid for inertial frame of 
reference only. The Euler angles 6,<j},y/ define the transformations that transform earth 
fix axes to body axes at any particular instant of time. The transformation is obtained by 
rotation of body axes, the order of rotation in which is very important. 



'cosy/ 

-sin^ 

0^ 

r cos# 

0 

sin#'' 

1 

0 

0 > 

V = 

earth fixed 

sin^ 

sin^ 

0 

0 

1 

0 

0 

cos^ 

-sin^ 


l o 

0 

K 

^-sin# 

- 

cos#. 

1 ° 

sin^ 

COS0 j 


body fixed 


(3.30) 


<j> — p + #tan#sin^ + rtan#cos^ 

(3-31) 

#= gcos^-rsin^ 

(3.32) 

rcos^sec# + # sin^sec# 

(3.33) 

The equations used to transform the body-axis velocity (u, v, 

w) into earth fixed axis are: 

X=ucosy/ cos# + v (cos^sin#sin^-sin^ cos^)) 

+ w (cos^ sin # cos <j> + sin y/ sin <f) 

(3.34) 

T=wsin^ cos# + v (sin^ sin#sin^ + cos^ cos^) 

+ w (sin^ sin#cos^-cos^ sin^) 

(3.35) 

Z-us\n9 - v cos#sin^~ wcos#cos^ 

(3.36) 
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3.2.3 WIND MODELLING 


The wind components of velocity are also taken into account by appropriate wind 
modelling. 

u'= u-W x cosO cos y/ + W y cos 6 smy/ + W x s\n6 (3.37) 

v : =v — W x (cos ^ sin 0 sin ^ - sin ^ cos <f>) - W y (cos^ cos^ + sin^ sin#sin^) 


+W, cos# sin ^ 

w'=w—W x (cos y/ sin 0 cos (j ) + sin y/ sin <f>)—W y (sin y/ sin 8 cos (f > — cos y/ sin 
+ W z cos 0 cos <f> 


(3.38) 


(3.39) 


The wind relative angle of attack and sideslip angles are computed using 

a = tan -1 (3.40) 

u 

I V 

P = tan' 1 — (3.41) 


3.2.4 STANDARD ATMOSPHERE MODELLING 

The International Civil Aviation Organization (ICAO) atmosphere, which is the 
most commonly used by aero ballisticians today is considered in the trajectory modelling 
in order to estimate the density and temperature variations with altitude. The ICAO 
standard atmosphere gives values of temperature and pressure and density up to an 
altitude of 20 km. It assumes sea level temperature of 15 0 C and a pressure of 101.325 
mbars and density of 1.225 kg/m 3 . The variation of density and velocity of sound are 
modeled using the following. 
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p = 1.225(1 -0.000022557695 .Z) 4 25587 


(3.42) 


^ = 340.249( 


288.15 -0.0065.Z os 
288.15 


(3.43) 


3.2.5 AERODYNAMIC MODEL FOR SIX DEGREE OF FREEDOM MODEL 


C>: — C Na M + C NB S + C Nq . 2 V 

(3.44) 

C,=C^ + C^+C,,.|t 

(3-45) 

C.=C^ + C, B S+C^ 

(3.46) 

Cy = C„./> + C^Sa + ■ C„.&- + C„ + C^r 

(3.47) 

C, =C, p .p+ C lBa Sa+C„,.^- + C lr .^+ C Wr Sr 

(3.48) 

c, = C„./J + cJ-^+C^Sa + C v Ac, lr Sr 

(3.49) 

Once the coefficients C N ,C r ,C x ,C m ,C n ,C, are suitably modeled then the six-degree 

of freedom model equations serves as a fairly accurate tool for trajectory modeling. 

3.3 DEVELOPMENT OF CONTROL INPUT STRATEGY FOR 

MISSION REQUIREMENT 


The mission requires that the missile should have a leveled flight at height of 
approximately 500 meters. Thus development of control input strategy through the fins so 
as to obtain the desired trajectory was an important step in the trajectory modeling of the 
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missile. By control input strategy we mean, that for the desired flight path what should be 
the deflections given to the fins as a function of time. This strategy once developed can 
be used to pre-program the fin deflection. In order to get the required fin deflections, the 
trajectory model was run for different sets of deflections. The trajectory' was simulated 
with these deflections. The mission requirements were: 

1 .The missile should attain a maximum height of around 500 meters in shortest possible 
time. 

2.After attainment of maximum height it should maintain a leveled flight up to the 
terminal phase. 

By varying the time intervals and magnitude of fin deflections three sets of fin 
deflections were chosen .The angle of attack a , flight path angle y , height etc were 
monitored. Thus by careful study of all the parameters the strategy was developed for the 
control deflection during unguided phase. 

Also it was required that angle of attack fluctuations should be minimum and to 
achieve zero flight path angles while maintaining a level flight at maximum height. So 
these deflections were linearly varied w.r.t time and thus the deflections were trimmed so 
as to: 

1) Maintain zero y angle during the level flight. 

2) Reduce the a fluctuations. 


50 



3.4 RESULTS AND DISCUSSIONS 


The fin deflections were given as a function of time and are shown in Fig 3.4.1— 
Fig3.4.5 The five set of different fin deflections shows the systematic approach to freeze 
the fin deflection requirement. These sets are tabulated in Table 3.4.1 
Table 3.4.1 Control deflections for different sets 


S.No 

tl 

3 

t2 

^2 

t3 


Set I 

0<t <5 

-10 

5<t<S 

+10 

8 < r < 12.5 

OH 


0<t <6 

-8 


+8 

9<f <12.5 

m 



-9 

5<t<9 

+7 

9<r <12.5 

-4 

Egg 

0 < t < 6 

-8 

6<t< 9 

+8 

9<f <12.5 

-3 

Set V 

0 <t <6 

-7 

6<t< 9 

+8 

9<r <12.5 

-1 


It should be noted that the fin deflections were given constant for Set I-Set IH and 
linearly varied between the time intervals tl-t2 and t2-t3 in Set IV to Set V. 

From Fig 3.4.1 it can be seen that a height to 500 m is not achieved. Also there are 
undesirable fluctuations in angle of attack, a . Thus another set of deflections was 
chosen so that it can reduce the maximum height attained by the missile. As we know 
that a negative fin deflection results in a pitch up moment so it was necessary to give a 
lesser fin deflection in tl interval in Set I. Fig 3.4.2 indicates that the missile is able to 
follow the desired trajectory path but it’s angle of attack fluctuation, is certainly not 
permissible. This pattern is also followed in Fig 3.4.3.The angle of attack, a fluctuations 
in all the above cases was due to the sudden change in deflections from negative to 
positive in zero time intervals. So it was necessary to apply the deflections gradually. 
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In Set IV & Set V fin deflections were applied gradually maintaining a linear variation 
from tl to t2 .Fig 3.4.4 shows the trajectory parameters obtained as a result of Set IV 
deflections. It is clear that the angle of attack, a got reduced substantially however 
requirement of 500 m height is not achieved. Further it was necessary to reduce fin 
deflections (<$,). 

Set V gave quite a promising configuration as shown in Fig 3.4.5. These set of 
deflections were frozen during all coarse of the subsequent analysis in the trajectory 
modeling. The final pre-programmed fin deflection used for further analysis is 

8 X = -7 deg. 

S 2 = +8 deg. 

<£, = - 1 deg. 

It can be seen that with this pre-programmed fin deflection, the mission requirement of 
reaching 500 m as well as having level flight at 500 m with reduced angle of attack 
oscillation is achieved. 
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th Set II deflections. 













Fig 3.4.4 Time history of flight parameters with Set IV deflections 




Chapter 4 


TERMINAL GUIDANCE AND CONTROL REQUIREMENTS 
4.1 GENERAL 

Theoretically the proportional navigation law issues acceleration commands, 
perpendicular to the instantaneous missile target line of sight, which are proportional to 
the line of sight rate and closing velocity 14 . Mathematically, the guidance law can be 
stated as 

a c = N'V C .X 

where a c is the acceleration command, N' a unit less designer chosen gain (usually 3-5) 
known as the effective navigation ratio, V c the missile target closing velocity and /l the 
line of sight angle (in radians). The over dot represents the time derivative of the line of 
sight angle or line of sight rate. 

In tactical missiles, proportional navigation guidance commands are usually 
implemented by moving fins or other control surfaces to obtain the required acceleration. 

4.2 SIMULATION OF PROPORTIONAL NAVIGATION IN TWO DIMENSION 

As per the law of proportional navigation, the normal acceleration required to 

engage the target is given by 

a c = N'.V c A (4 ' 1} 

X is computed using the following relation 

■ rxV (4.2) 

f/V — ~ ’ 
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where 


r =(x T -x).i +y.j 


V — {x T — x).i + y.j 


X can be expressed explicitly as 


V c , the closing velocity is computed as 


(4.3) 




[(x r -x)(-i) + y] 


(4.4) 


c ((x r —x) 2 +y 2 ) m 

The normal acceleration was then resolved along inertial x and y directions and 
solved for x, y to complete engagement simulation using following mathematical model. 


d 2 x 

dt 2 

£ 

dt 2 


= -a n sin X 


= a n cos X 


(4.5) 

(4.6) 


where 


— = V r =Vcosy 
dt x 

— = V v -Vsiny 
dt y 

where y is the flight path angle and a c is computed using eq (4..1) 
Initial conditions used for the simulations are 


(4.7) 

(4.8) 
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At t — 0 , X — Xq y — 

X = (V COS^) /0 
j> = (Fsin 0) t=o 

These equations were solved using fourth order runge kutta routine to arrive at missile 
trajectory 

4.3 GENERATION OF REQUIRED NORMAL ACCELERATION 

Based on kinematics analysis, it was found that to engage the target a definite 
value of normal acceleration (a c ), perpendicular to the line of sight, must be generated in 

flight. As stated earlier, this acceleration needs to be generated by the aerodynamic 
normal force acting on the missile. To ensure that sufficient acceleration is produced, a 

# 

definite value of angle of attack must be introduced to the missile. The angle of attack is 
introduced by moving the tail in particular direction. The missile doesn’t instantaneously 
respond to tail deflection as far as generating steady state acceleration is concerned. It 
takes definite time to build the trim angle of attack corresponding to a definite value of 
tail deflection. Since the missile is having low damping, it is must to device a stability 
augmentation system to artificially damp the oscillation so that steady state angle of 
attack, thus the acceleration is obtained almost instantaneously. 

In this section, study has been carried out to model the dynamics including the 
expected value of trim angle of attack using non-linear aerodynamic model. Finally a 
tentative stability augmentation system has been proposed to enhance damping to obtain 
fast steady state response of the missile. 
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4.4 ESTIMATATION OF TRIM ANGLE OF ATTACK AND NORMAL 
ACCELERATION USING NON LINEAR MODEL 


At trim angle of attack, net moment on the missile must be zero. The expression 
for moment coefficient for body fin with tail deflection is given by 

^ m bodyjm m ) body ,nonhn 


hail 


(C) 


bodyjm 


l.{k 2 -k x ) + rjCd c .-^-. sin 2 a 


A 




-Hep, body) 


(4.9) 


(4.10 ) 

Following approximation was made to arrive at approximate analytic expression for trim 
angle of attack. 

1. sina . 

2. Centre of pressure of body and fin were assumed to be independent of angle of attack. 
The moment coefficient C m can be expressed as 


C m =Y v a + Y 2 .a 2 +Y r S 


^1 (^'Na bodyjm 


{■^■Cg X cp.hody ) , (r ^ {-^cg X cp j „ ) 

• ; j 


Y 2 =rjCd c . 


■^■p (*- ^ cp, body ) 

A-nf d 


y -(r \ fe X cp j m ) 

Z 3 ~V^Na'fm' jr 


(4.11) 

(4.12) 

(4.13) 


At trim ,moment coefficient is zero, therefore, for a given fin deflection, 5 , we get the 
equation for the trim angle of attack a mm to be 
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0 -^ - a trm + *2 - a tnm +YyS (4 J 5) 

we can use quadratic formula to solve for the trim angle of attack. After eliminating the 
unrealistic root we get 



(4.16) 


The expression a tnm for different tail setting angle, 8 has been used to compute steady 
state normal acceleration .It may be noted that at |a B | = |a c | if 0 = 0 ,i.e. when missile 
aligns with the instantaneous line of sight. 



(4.17) 


4.5 DESIGN OF STABILITY AUGMENTATION SYSTEM 

Schematically the stability Augmentation system may be represented as follows: 



Fig 4.5.1 Schematic of a Stability Augmentation System (SAS) 

The rate gyro flight control system artificially increases the low damping of the 
open loop flight control system by the use of rate gyro sensor and the principle of 
feedback. The tail deflection is proportional to the pitch rate to enhance pitch damping. 




62 





Pitch damping ratio, £ and natural frequency <x> n obtained by pure pitching 


approximation 11 is: 


II 

N> 

in i 

(4.18) 


(4.19) 


For short period approximation : 


-(^,+f-) 


2 JM a .=2—M a 

V q v 


co n =AM q .-^--M a 


where 


1 rrl 


M, 


_ 2 


pV^SdC. 


yy 




{pV'SdC,, , d , 
\2V, 


~(Cu+cJ± P r 2 s 

Z„ = 2 

m 

Time to half amplitude and angle of attack response is computed using 


l V2 


0.693 

%.co n 


a(t)=a tr J+-j==^sin(^h-fcy n t + £) 


(4.20) 


(4.21) 


(4.22) 

(4.23) 
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(4.24) 


a. 


trim 


K 


e = tan 



(4.25) 


For a pre-decided T y2 amplitude, £ required is computed using expression for 
short period approximation as given in eq (4.22) & (4.20) respectively . In estimating <»„ 
using eq(4.21) , M q is assigned zero value. Once £ is known, for a given value of 
Z a ,M a ,V the required M q value is estimated using eq (4.21). The derived M q is 


attained artificially using ‘ q ’ feedback to tail ( A S e =kq). 

Change in moment coefficient, A C m , due to additional tail deflection is calculated using 


— C mSe AS e 
Since A S e = kq 


Thus the gain, k is estimated as 


k = 



required 



'mSe 


d 

2V 


(4.27) 


The tail when deflected by kq , results in reduction inT 1/2 amplitude for a design value of 


£ (required). 
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4.6 RESULTS AND DISCUSSION 

4.6.1 TWO DIMENSIONAL TARGET ENGAGEMENT SIMUMATION- 
NORMAL ACCELERATION REQUIREMENT 

For target engagement simulation, the target was assumed to be at 2000m.The 
missile was activated for guidance at a height of 500m as geometrically shown 
in Fig 4.6.1. Missile was fired at various angles ranging from 0 = - 20° to +30° at a 
velocity 290 m/sec. Fig 4.6.2 pictorially presents target engagement simulation. 



Fig 4.5.1 Two Dimensional Missile Target Engagement geometry 

Target engagement simulation was carried at for various values of (3-5). It was 
found that for N' >= 3, the simulation results in exact matching. Though literature 
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recommends value of N' to be with in range of 3-5, however from the point of view of 
practical limitation, it is desirable to keep N'~ 3. 

Results were generated for different angle of projection towards the target. The 
value of iV'=3, ensures exact engagement of target for different angle of projection 
( 0 = —20 deg to +30 deg). This is pictorially presented in Fig. 4.6.2.however acceleration 
requirement for each case will be different for different values of initial launch angle, 
0 . It is obvious that for 6 = 30 deg ,it will require largest normal acceleration to engage 
the target. In order to estimate the magnitude of normal acceleration (a c ), based on 
proportional navigation kinematics, the code was run at V=290 m/sec for different values 
of 0 and a c was recorded as presented in Fig. 4.6.3. It can also be seen that for 9— 0, i.e. 

the missile fired horizontally, the a c required is around 3g and for #--10 deg (relatively 
closer to the line of sight) the a c required is around lg.In practice ,since the missile is 
terminally guided , 6 should be set close to line of sight with small leading angle, thus 
acceleration requirement can be reduced accordingly. 
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Height ofRelease.m 



range, m 


Fig 4.6.2 Two-dimensional missile target simulation 
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Acceleratlonfm/s 2 ) Normal to L.O.S 



Fig 4.6.3 Variation of required acceleration normal to LOS for different 
release angle 6 . 
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4.6.2 ESTIMATION OF TRIM ANGLE OF ATTACK AND NORMAL 
ACCELERATION 


Trim angle of attack and normal acceleration are computed using eq (4.16) & 
(4.17). The values of trim angle of attack and normal acceleration a n are computed for 

different tail setting angle and velocities. Table 4.6.1 presents the estimated values of trim 
angle of attack and normal acceleration for different setting angle and velocities. Figure 
4.6.4 pictorially presents the result. 

Table 4.6.1 Theoretical Estimation of a mm & Normal acceleration a n for 


different velocity and S input. 


S , Deg. 

Velocity 

<**m 


5 

200 

-5.316 

4.3294 

5 

230 

-5.316 

5.7256 

5 

260 

-5.316 

7.3166 1 

5 

280 

-5.316 

8.4855 

5 

290 

-5.316 

9.1025 





10 

200 

-10.201 

10.529 

10 

230 

-10.201 

13.925 

10 

260 

-10.201 

17.794 

10 

280 

-10.201 

20.637 

10 

290 

-10.201 

22.137 





20 

200 

-19.012 

26.733 

20 

230 

-19.012 

35.355 

20 

260 

-19.012 

45.179 

20 

280 

-19.012 

52.397 

20 

290 

-19.012 

56.207 





' 25 

200 

-23.047 

36.076 

25 

230 

-23.047 

47.711 

25 

K 260 

-23.047 

60.969 

25 

r ~~ 280 

-23.047 

70 709 

25 

^ 290 

-23.047 

75.85 





30 

200 

-26.883 

45 912 

30 

230 

1 -26.883 

60 719 

30 

' 260 

-26.883 

77.592 

30 

280 

-26.883 

89.988 

30 

290 1 -26.883 

96.531 
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Velocity, m/sec 

Fig 4.6.4 Variation of Theoretical a n at a inm for different fin deflection 

Based on this study, it can be stated that to generate a 3g acceleration, the tail deflection 
of around 1 5 deg is required at a velocity of 290 m/sec. 


4.63 DESIGN OF STABILITY AUGMENTATION SYSTEM (S.A.S) 

To estimate the value of gain ,k for updating pitch damping using pitch rate 

feedback the eq (4.27) has been used. (C mq ) ^ was calculated for a flight condition of 

velocity=290 m/sec at a height of 500 m. Table 4.6.2 lists the values of | required for 

different desired Ty 2 amplitude. 
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Table 4.6.2 Values of T m at Vel=290 m/sec and Height=500m 


^1/2 ( Sec ) 

£ (Required) 

k 

0.1 

0.427 

0.04 

0.05 

0.854 

0.08 



Fig 4.6.5 Variation of a damping for different SAS coefficient values. 


For a specific case, as stated earlier 7^ /2 =0.1 sec is considered and the feedback 

gain k was estimated. Now six degree of freedom model was used to compute angle of 
attack time histories with and without feedback. For this specific case, the angle of attack 
response without feedback assumed zero pitch damping in open loop. 


71 





Fig 4.6.6 A comparative time history of a with and without feedback 

A comparative time history of a with and without feedback is presented in Fig 4.6.6 .It 
can be easily verified that the steady state value reaches almost instantaneously. Thus it is 
validated that the scheme for stability augmentation system will enhance pitch damping 

to designed value. 



Fig 4.6.7 variation of fin deflection for different SAS coefficient values 
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CHAPTER 5 


CONCLUSION AND SCOPE FOR FUTURE WORK 

5.1 CONCLUSION 

1. The task of generating aerodynamic parameters using theoretical wind tunnel methods 
have been completed. Based on the comparison it is observed that the functional form 
used in theoretical methods can be used to generate initial design parameters. The 
alteration of flow field due to shedding vortices could not be captured using theoretical 
methods, thus for fine tunning, wind tunnel estimates should be used. 

2. Trajectory model has been developed to pre-program the tail deflection to satisfy the 
requirement of reaching height of 500 m and maintaining level flight up to 3000 m. 

3. Normal acceleration (perpendicular to the line of sight) required to simulate target 
engagement has been modeled based on proportional navigation kinematics. 

4. Missile response to generate the desired acceleration per unit tail deflection has been 
evaluated. Model to estimate expected steady state angle of attack (a tnm ) has been 
proposed. In order to reduce the T m amplitude of missile airframe response, stability 
augmentation system has been proposed. 
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5.2 SCOPE FOR FUTURE WORK 


1 .Exhaustive wind tunnel testing for angle of attack sweep to fixed tail deflection should 
be carried out for fine-tuning. 

2. Theoretical methods need to be updated to capture ‘effect of shedding vortices’ and fixe 
lifting characteristics of the missile. 

3. Target engagement simulation using proportional navigation needs to be developed 
using six degree of freedom modeling. 


74 


REFERENCES 


1 . Kennedy , G.P. /‘Missiles and Spacecraft of National Air and Space 
Museum”, Smithsonian Institution Press Washington DC, 1983. 

2. Benecke,T.,and Quick, A. W. ,” History of German Guided Missile Development,” 
Proceedings of AGARD first Guided Missile seminar, 1956. 

3. Nesline, F.W. , and Zurchan , P. , “A New Look at Classical Versus Modem Homing 
Guidance, ’’Journal of Guidance and Control, Vol.4,Jan-Feb,1981, pp.78-85. 

4. Yuan , C.L. , “Homing and Navigation courses of Automatic Target-seeking 
Devices,” RCA Labs.,Rept. PTR-12 C , Princeton , NJ,Dec 1942. 

5. Benett , R.R. , and Mathews , W.E. , “Analytical Determination of Miss Distance for 
Linear Homing Navigation Ststems ,” Huges Aircraft Co. ,TN-260 , Culver city , CA„ 
March 1952. 

/ 

6. Fossier , M.W. , “The Development of Radar Homing Missiles,” Journal of Guidance 
Control , and Dynamics, Vol. 7 ,Nov-Dec. 1984, pp. 641-651. 


75 


7. Yuan, C.L., “Homing and Navigation Courses of Automatic Target-seeking Devices,” 
Journal of Applied Physics, Vol ,19,Dec.l948 ,pp. 1122-1128. 

8. Bryson, A.E., and Ho, Y.C., Applied optimal Control, 13 Laisdell, Waltham, MA, 
1969. 

9. Joset Rom , High angle of attack Aerodynamics , Subsonic ,Transonic and Super sonic 
Flows ,springer-verlag New York , Inc, 1992. 

1 0. Moore , Frank .G. /‘Approximate Methods for Weapon Aerodynamics”, Progress in 
Astronautics and Aeronautics, Apr 2000. 

11. Maine, R.E., and Illiff, K.W. , “Application of Parameter Estimation to Aircraft 
Stability and Control,” NASA RP 1 168,1986. 


1 2. The Engineering Design Handbook, “American Material Command, AMCP,706 -280. 


13. Hoak , D.E. , USAF Stability and Control DATCOM; Air Force Flight Dynamics 
Laboratory, Wright Patterson Air force Base, Ohio, 1960, revised 1975. 


14. Zarchan, P., “Tactical and Strategic Missile Guidance,” 3 rd AIAA, Inc.1947. 


76 


APPENDIX ‘A’ 


data BASE FOR ESTIMATION OF (C„)„„ , (C,)^ etc. 

Tabic 2.2.1. (a) (C„) MJ , as a function of Mach number, aspect ratio & taper ratio 
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DATA BASE FOR ESTIMATION OF INTERFERENCE FACTOR 

TABLE 2(B) 
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Data for a D at <p = 0 deg 
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APPENDIX ‘B’ 

PROCESSING OF WIND TUNNEL DATA 

Two sets of test were undertaken 
1 .Dry run 
2. Wind run. 

Both these sets were carried out for three configurations: 

1 .Body alone 

2. Body + fin for delta=0 

3. Body 4 fin for delta=-40,-35,-30,-25,-20,-10,+10. 

The dry run for all the above three configurations are considered offsets and are 
subtracted from each of the values of wind run. 

Thus a new set of corrected values of cl,c2,c3,c4,c5,c6 in voltage form is obtained 
which is processed to obtain forces and moments and their respective coefficients. A 
detailed procedure is described below. 

The wind tunnel data was obtained at the balance center in the form of 
cl,c2,c3,c4,c5,c‘6 in voltage form. Also the slope was obtained. 

If 

y-output from wind tunnel in voltage form. 
m=slope(given) 

x^force (kgf) or moment (kgf-m) 

Thus the values of forces and moments axf , n\, n2, s\,s2,rmvn. kgf and kgf-m respectively 
were obtained by solving the equation : 
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y=mx 


Thus following data was obtained w.r.t the balance center. 
Sign Convention(For tunnel axes) 

1 ) positive x is towards tail . 

2) Positive y is rightward as viewed from rear. 

3) Positive z is upward. 

Sign Convention(For Trajectory axes) 

1) Positive x is towards Nose. 

2) Positive y is rightward as viewed from rear. 

3) Positive z is Downward. 

Also it should be noted: 

f x =-f a 

Step 1. 

Force (balance center) 

i^x) balance = ~ CtX f 

(K)u^, = < nl+n2 1 

W )«-. = *!+ *2 

Moment (balance center) 

(^w=-(' !l - n2 >* 0065 

(Y m )/Hhmx =(sl-s2)x 0.065 
(*»)««. = rm 
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Step 2. 

It should be noted that the above values are w.r.t the balance center. These needed to be 
transformed to the model center of gravity. 

The distance between the balance center and the model eg was calculated and thus if xd 
is the distance of balance center and model eg then the following equations are used to 
transform the force and moments from balance center to model c.g 
Force Transformation from balance center to model center of gravity 


balance 




balance 




balance 


Moment Transformation 


from balance center to model center of gravity 


( J>m )c K = ( Pm ) Man c< ~{F N ) ha[ance .xd 


(Ym\,, =(Ym) Monce ~{F r ) halance xd 


=(Rm) lHllimct 

Step 3. 

The next step was to evaluate the forces and moment coefficients. 

(F t ),,.9.Sl 
' V ” 0.5 pV%., 

c (^)<r 9 - 81 
T 0 .5pV 2 S K/ 

c _ 

* Q.5pV%, 
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c (-Pm) g .9.81 

- 0 .spvXf-J 
C = J^k!£L 

” 0.5 pV'S^d 

c W.,-9.81 

Step 4 

The next step is to transform the body axes coefficients to wind axes . 

Force Transformation from Body to Wind Axes 

D = {F X ) cg cos a cosy/ + (F N ) cg sin a cos y/ + (F y ) cg sin yr 

/, = ( F n ) cg cos a cos ^ - (F x ) cg sin a cos <f> - (F Y ) cg cos yr sin tf> + (F x ) cg sin yr sin <f> 
Y = ( F y ) CK cos yr cos <j) - (F x ) cg sin y/ cos - (F w ) cg cos a sin <j> - (F x ) cg sin a sin ^ 

Moment Transformation from Body to Wind Axes 

Pm = ( Pm) cg cos y/ cos <j> + (Rm) cg sin y/ cos a cos tj> + (Ym) cg sin a sin ^ cos <j) 
-(Ym) cg cos asin<f> + (Rm) cg sin a sin (f> 

Km = ( Fwj) ( . s cos cccos<f>- (Rm) cg sin a cos y/ + ( Pm) cg cos yr sin <j> 

+(Rm) cg cos a sin (f>smy/ + {Ym)^ g sinasin^sin^ 

Rm = (Rm) cg cos a cosy/ + (Ym) cg sin a cos yr - (Pm) cg sin yr 

Step 5 

Conversion of Wind Axes Forces and Moments to Coefficient Form: 
Force Coefficient 

Zx9.81 
= O.SpV 2 S ret 
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Dx 9.81 


' D 0.5 pV 2 S ref 
Kx9,81 
' K " 0.5 pV 2 S rcJ 

Moment Coefficient 

_ Pmx.9.%\ 

Lmx9.81 
' n ~ Q. 5 pV 2 S ref .d 
i?mx9.81 
'' [).5pV 2 S ref .d 


Thus body axes coefficients C x ,C Y ,C N ,C m ,C n ,C l and wind axes coefficients 
C p , , ( 7 , , C ' m , C n , C, are obtained for each of the configurations discussed previously as 

a function of angle of attack. 

Let us now denote the body axes coefficients obtained after processing the raw wind 
tunnel data for three configurations namely: 

1. body alone 


^ " N )aixA' 
X ) hotly 

)htxlv 

KvU, 



2. (body + fin) M 

N \h,xly ~ ^ N 
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(CA 


c v 


v\hody-{ fm) 6 4 ~~ Cy 

(C ) =zC 

v "7 n body + fin ) 6A ^ni 

n \hodyk fm) A () ~~ 

^ / \bodyi fm) d 0 


3. (body + fin) & 


s^s\ 


N\hodv¥jm\ USv C N 

^ X ) { body 4 fin ) A Al ~~ C x 
(^ Y )( Wy4 fm , ~ C y 
(C m \ body + fin ), v <VI ('m 
(^ « ){ body* fm ) 6 Al ~ ^'n 
^ / \body \ fm ) A Al C i 


Flic next step was to develop a mathematical model, which can evaluate the stability 

derivatives 

Longitudinal derivatives : C Na > C N $ , C N(J , C ma , C w<? , C m<5 
Lateral Derivatives: C y p , , C yp , C^ r , C n p , C wr , C w;7 , C n ^ r 

Cip>C lSa ,C lp yC lr 

Since the test was carried out for 
1 .Body alone 

2. Body 4* fin for delta=0 

3. Rody + fin for delta=-40 ? -35,-30 r 25 r 20 r 10,+10. 

Thus we have assumed that all the stability derivatives are composed of the contribution 

from: 

i )hody alone 
\\)(body + fin) Sm 0 
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iii) ( body + fi n )s*s\ 


Thus 

c N =(C,) w +(C N W +CQ/W. 
c x = (C x ) body + (CxW + W« 

C y = (C y )** + ( c r)(/«),=o + 

c m = ( c m ) body + (C.W + (C - W. 
c, = (C,U +( c 'X^w +( c »>o«»)«i 
c„ = (c„) body +(C„\ fi ^ o +(C « w. 

Stability Derivatives From Body Alone ^ &om the 

In order to evaluate (C# a )&**,,>( rp)body 

farther resolved mto individual components. 

wind tvmnel in coefficient form was further resoi 

(C,/),,,,,, =( c «»)*»*’ a mnment coefficients 

„ the body axes force coefficients and moment 
Since we have the bony a* 

r r C C C ,C, from the wind tunnel. 

(fi Nobody 

(pNa^body ~~ ^ 

(Pm) body 


a 


( Pma Podv 

cc^w-^u 
(C,„) h ,,» = -2-( c »«W - T 

Stability derivatives from Body + fin at <5 = 0 

fC ) ^(CP)(body+fM)s*4 m b ° dy 


88 



Let (C N ) (/in)s ^ - {C N ) fm 


(C Na ) 

(C ma ) 


(C N ) 


fin 


'fm 


a 

(CJ 


'ma J fm 


m) fm 


a 




=-2 -(C m ) rm . 


ma / fm 


( xcg-xcp ) 

fin 

d 

f xcg-xcp 

fin 

V * " 


(Cyp) ^ = ~(C Na ) yj„ 


(^). = ~2-(C~) . 


< C A — 2.(C*) V 

C n ,=-2.15.(C„ a ) /01 . 


(xcg-xcp ) 


i? 


'r 2 


Stability Derivatives From Body + Fin at 8 = 81 (-40, -35, -30, -25, -20, -10, +10) 

= ^N){body + fm) Ml ~(C N )bcdy 

~(C m )(kody+fin) SmSl — (Qn) (/lnfe) — (Qi 

Let 

81 = 8 

= (C y ), 

=(£>,)(? 

r _ 
w £1 

_ (Q»)<? 

'-'mS ~ ~ 
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c =-r 

^ \>Sr v NS 


C = -C 


r* =C 

^ nbr mb 

In sum the total derivatives discussed can be written in summation form as 

^ Na ~~ <V«r ) /?» 

^ ma ma^bodv //« 

^ mq “ wt/ )/>w/v (^wc/ ) /w 

^ ' M/ A/<r/ )a«/v % ) fm 

( vft \'fi)hodv (C vf ) f m 

^ v/ = vr ) fm 

^ vp = (X' yp )bodv ) fm 

^ = (C'nfi )body ^nfi) fm 

Futher 


r =C 


Q- = 2.C V/J . 


' xcg-xcp 


\ 


^ ' / f >V War ) /w ’ 


Ida ~~ iV« ) fm ' 


'IY 

£ 

d 


Thus the final model used in the trajectory modeling including the damping derivatives 


— ^ A/o'® + + 2 V 

^',4 = ^/i« - a + ^ + 2v 

(7(5? 

Cm = £*»!« - a + ^nxS + ^mq ' *2^ 

»(5? 

c, = C v p.fi + C vS .S+C >p .-^ 


+ C yr 


rd 

2v 


+C 


d> 


2v 
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C^C^+C^c+C^+C^+C^r 

C.=C^.fi+C^+C^+C^+C^ 

These derivatives can be used in the trajectory modeling. The derivatives are obtained as 
a function of angle of attack. 
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